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Abstract 

High precision multipole polarizabilities, a£ for £ < 4 of the Is ground state of the hydrogen 
isoelectronic series are obtained from the Dirac equation using the B-spline method with Notre 
Dame boundary conditions. Compact analytic expressions for the polarizabilities as a function 
of Z with a relative accuracy of 1CP 6 up to Z = 100 are determined by fitting to the calculated 
polarizabilities. The oscillator strengths satisfy the sum rules £^ f$ = for all multipoles from 
£ = 1 to £ = 4. The dispersion coefficients for the long-range H-H and H-He + interactions are 
given. 

PACS numbers: 31.15.ap, 31.15.ac, 34.20.DK 
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I. INTRODUCTION 



The present paper reports calculations of the polarizabilities of the hydrogen atom and 
isoelectronic ions using the Dirac equation to describe the underlying dynamics. Such calcu- 
lations are now topical since some atomic polarizabilities can directly impact the definitions 
of two fundamental quantities, the Kelvin and the second The new generation of optical 
frequency standards have reached such precision that they are sensitive to the black-body 
radiation of the apparatus itself 0|. The resulting black-body radiation shift is largely de- 
termined by the differences in polarizabilities of the two atomic states involved in the clock 
transition. Additionally, very high precision measurements of the helium dielectric con- 
stant have been recently reported |3j. In conjunction with high precision calculations of the 
static dipole polarizability [4|, these measurements can result in improved determinations of 
Boltzmann's constant and thus the Kelvin 0, 0|. 

Another reason for doing such calculations is that they can be used to verify the accuracy 
of computational methods and tests of fundamental theory. The polarizabilities of the 
hydrogenic ions are properties of the ground state of a set of systems that are often used 
to test the fundamental principles of physics. It is rather surprising that the first explicit 
calculations of the quadrupole polarizabilities of the hydrogenic ions based on the Dirac 
equation have only just been reported [fj]. 

An important advance in the topic of the dipole and higher multipole polarizabilities 
was an investigation based on the Pauli approximation that gave expressions for the static 
multipole polarizabilities up to 0(a 2 Z 2 ) [7J. This was a generalization of an earlier work 



which gave the 0(a 2 Z 2 ) expression for the dipole polarizability |8|. The work on dipole 
polarizabilities was extended tofaZ) 4 {aZ) 6 [Io[ and to all orders in terms of a general- 



ized hypergeometric function 11] . Apart from a very recent calculation [6|, the expressions 



for the quadrupole and higher-order polarizabilities have not had independent confirmation. 
There have been a number of independent calculations of the dipole polarizability and re- 
lated sum rules. Many of these investigations have been computational in nature. Drake 
and Goldman derived expressions for some dipole oscillator strength sum rules as well as 
performing some explicit calculations of the dipole polarizability [12| by ex pan ding the wave 
function as a linear combination of exponential type functions. Goldman [l3( extended the 
basis set approach to calculate the dipole polarizability of hydrogenic ions from Z = 1 to 
Z = 115 using a Gauge-invariance method. A fit to the calculated polarizabilities was used 
to create an (aZ) n expansion of the polarizability including terms up to (aZ) 8 . There 
have been a number of other computational investigations of the dipole polarizabilities of 
hydrogenic ions based on Dirac equation (I3-20l . 

The present calculations used the B-spline Galerkin method with Notre Dame (ND) 
boundary conditions 19 1. Other approaches to the B-spline boundary conditions have been 



proposed [l7|, |l8|, |2l|-|23|. There is at present no overwhelming reason for adopting more 



complicated boundary conditions in preference to the ND boundary conditions. The B 



spline approach to atomic structure has a number of advantages |25| , it does not lead 



to linear dependence, the basis can be made effectively complete in a finite region of space, 
the details of the basis are easily adjustable and results are numerically stable. However, 
like all basis set approaches, the energy spectrum also has a sea of negative-energy states 
(the Dirac sea) and it is also possible for spurious states to appea r in the positive energy 
spectrum 18j]. These issues have been discussed extensively 



The present B-spline calculations of the multipole polarizabilities give numerical values 
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that are more precise than any previous calculation. Values of associated oscillator strength 
sum-rules are also given. The nuclear mass was set to be infinite and the point nucleus 
model was adopted. Values are reported for intermediate sums including the entire set of 
states and also for a set of calculations that omitted the negative energy states from the 
Dirac Sea. Analytic expressions for the polarizabilities are constructed that are accurate to a 
relative precision of 10 -6 for Z < 100. The static multipole polarizabilities for quadrupole, 
octupole and hexadecupole transitions have been computed and found to be compatible 
with the 0(a 2 Z 2 ) expressions of Kaneko et al [7]. The sum rules Ylifof = 0' P rov ide 
a valuable consistency check on the reliability of our calculations. Finally, the dispersion 
coefficients that describe the long-range interaction of the H-H and H-He + dimers in their 
ground states are presented. All results are reported in atomic units and the value of fine 
structure constant, 1/at = c = 137.035 999 074 [26| was used in all calculations reported in 
this work unless specifically mentioned. 



II. FORMULATION 

A. Dirac equation of single-electron atomic system 

The single-electron Dirac equation is 

H D ^(r) = E^{v) , (1) 

Hp is the Dirac Hamiltonian, 

H D = ca-p + I3mc 2 + V(r) , (2) 

where m is the electron mass, c is the light velocity, p is the momentum operator, cx and /3 
are 4x4 matrices of the Dirac operators 0]. 

The wavefunction for the hydrogen-like ion can be written 

■0(r) = - ( ^ nK ( r )^ Km (^") j (3) 



r 



where P ntt (r) and Q nK (r) present the larger and small components of radial wavefunction, and 
Q Km (r) and f2_ Km (r) are corresponding to the angular components. The angular quantum 
number k are connected with j and £, 

n = £(£ + l)-j(j + l)-l/A, (4) 

Substituting Eqs. (T5]) and ([3]) into Eq. ([1]) and separating the radial and angular components, 
gives the following coupled first-order differential equations for radial components P nK {r) and 

Qnn ('") j 

[V(r) + mc 2 }P nK (r) +c (J--^\ Q nK (r) = EP nK (r) , (5) 
-c + ^) P nK (r) + [V(r) - mc 2 ]Q nK (r) = EQ nK (r ) , (6) 
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In this equation V(r) is the interaction potential between the electron and nucleus, 



V(r) = ~, (7) 

with Z being the number of nuclear charges. 

In order to compare with non-relativistic calculations, we replace the energy E by e = 
E — mc 2 , and the radial Dirac equation can be written as matrix style, 



-c(f r + f ) -2mc 2 + V(r) J I Q nK {r) ) I Q nK (r) 



(8) 



B. B-spline Galerkin method 

The radial wavefunctions P nK (r) and Q nK {r) are expanded in a iV- dimensional basis of 
B-splines of order k, 

N 



P(r) = 5>B*(r), (9) 

i=l 
N 

Q(r) = ^^(r), (10) 



where the subscripts n, k have been omitted from the functions P„ K (r) and Qn K (^) for no- 
tational simplicity. The function Bf(r) only take nonzero values for the knot intervals 
ti < r < ti+k- The normalization condition is 

POO 

/ [P 2 (r) + Q 2 (r)]dr = l. (11) 
Jo 

The details of the B-splines and ND boundary conditions have been discussed in detail 
elsewhere 18, The large and small radial components are independently expanded in a 
B-spline basis with the boundary conditions, P(R) = Q(R) and P(0) = 0, where R is the 
radius of confining cavity. 

B-splines of k = 9 order were used with the endpoints of multiplicity 9. An exponential 
knot distribution for the B-splines is adopted, e.g. 

exp( 7 (^i T )) - 1 
U+k-i = R x ^""i-i" (12) 
exp( 7 ) - 1 

where i — 1, 2, • • •, Nx and N± = N — k + 2 being the maximal value of i. The exponential 
knot parameter 7 depends on the radius of confining cavity R, 

j = G(Z)xR. (13) 

The function G(Z) for Z > 2 satisfies the recurrence relation 

G(Z) = G(Z-1) + ^, (14) 
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where G(l) = 0.055 is an optimized value for the hydrogen atom. The confining cavity 
radius R (which is different for different Z) was chosen to reproduce the exact ground-state 
energy [27| of the hydrogen-like ions to at least 20 significant digits 

Exact 2 



where n is the main quantum number. 



(aZf 



n 



\K\ 



(aZ) 



212 



— C 



(15) 



C. Polarizabilities for the single-electron atoms 



In an weak external electric field, the static 2 £ -pole polarizability for an atom is usually 
defined in terms of a sum over all intermediate states including the continuum, 



At) 

J gi 



(16) 



The initial state, ip g {r), with energy, E g , is excluded from the summation over i, The 2^-pole 
oscillator strength f^f from ground state g to excited state % is defined 



(2*+l)(2j, + l) ' 

where j g is the total angular momentum for the ground-state. The wavefunction and energy 
of the excited states are ipi(r) and E%. C^(r) is the £-order spherical tensor. 

Using Eq. ()3]), the radial and angular parts of matrix element in the Eq. (ITT]) are 



(^ s (r)|r%(r)) 



(18) 



(0 Kff (r)||CW||^(f)) = (2^ + 15(2^ + 1) 



3g Ji t 

-1/2 1/2 



(19) 



Polarizabilities that are computed including both the physical states and negative energy 
states of the Dirac sea in Eq. f)16p are denoted by af. Polarizabilities that are computed 
by omitting the negative energy states of the Dirac Sea in Eq. (1TB"]) are denoted by a[ . The 
states of the Dirac sea are energetically distinct from the physical states The polarizabilities 
computed using the 0(a 2 Z 2 ) expressions of 0] are denoted as af. 

The polarizabilities can be expanded as a series in powers of (aZ) 2 . The series is written 



= a™ 



l + J2 X ^(aZ) 



2i 



i=l 



(20) 



where the non-relativistic multipole polarizabilities, af R , for the ground-state hydrogen-like 
ions, which have the exact values 128] 



a 



NR 



(2l + 2)!(l + 2) 
£(£+ l)2^Z 2i + 2 ' 



(21) 
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D. Oscillator strength sum rules 

There are a number of oscillator strength sum rules besides those which define the mul- 
tipole polarizabilities. We make the definition 

#(*) = £/j?(^ ( 22 ) 

i 

The expression with t = 1 and n = — 2 is the dipole polar izability, The case when I = 1 and 
n = is called the Thomas- Reiche-Kuhn (TRK) sum rule. In the non-relativistic calculation, 
Si(0) should be equal to the number of the electrons. The case with £ = 1 and n = — 3 is 
related to the non-adiabatic dipole polarizability j29[. One finds that S\(— 3) = 43/ (4Z 6 ) for 
non-relativistic hydrogenic atoms. The Si(— 1) coefficient [3(| is related to the long-range 
atom wall dispersion coefficient j3l[. One finds that Si(— 1) = 2/Z 2 for non-relativistic 
hydrogenic atoms. The relativistic sum rules are useful in testing the completeness of basis 
sets for variational representations of the Dirac spectrum j32j and set a foundation for testing 
other methods. 

As with the polarizabilities, the sum-rules can be evaluated by summing over all states, 
or just the positive energy states. Sum rules that are computed including both positive and 
negative energy states in Eq. (I2"2"j) are denoted by Sf(n). Sum rules that omit the states 
from the negative energy sea from the sum are denoted by S^{n). 

III. RESULTS AND DISCUSSIONS 

A. Polarizabilities and sum rules for hydrogen 
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FIG. 1: (color online). The convergence of ground-state energy (a.u.) relative to the exact Dirac 
equation energy of the hydrogen ground state. The number of B-spline basis functions is N, while 
the radius of confining cavity is R = 400 a.u. 

The difference of the B-spline ground-state energy from the exact energy given by Eq. ( fl5l) 
(this is -0.500 006 656 596 553 596 900 786 4298 a.u.) as a function of the dimension of the 
B-spline basis is plotted in FigJTJ This calculation was performed with a confinement radius 
of R = 400 a.u.. This ensures that none of the atomic sum rules reported in this paper are 
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TABLE I: The convergence of the static multipole polarizabilities (units of a.u.) for the hydro- 
gen atom ground state as the dimension, N, of the B-spline basis set was increased. The radius of 
the confining cavity is R = 400 a.u. 



N 










4 


4 


4 





100 4.499 751 495 18 14.998 829 821 131.2379 

150 4.499 751 495 177 64 14.998 829 822 856 41 131.237 821 447 83 2126.028 674 4992 

200 4.499 751 495 177 639 27 14.998 829 822 856 441 76 131.237 821 447 844 63 2126.028 674 499 1281 

250 4.499 751 495 177 639 267 48 14.998 829 822 856 441 70 131.237 821 447 844 661 2126.028 674 499 128 81 

300 4.499 751 495 177 639 267 398 14.998 829 822 856 441 699 67 131.237 821 447 844 662 144 2126.028 674 499 128 831 

350 4.499 751 495 177 639 267 396 1 14.998 829 822 856 441 699 61 131.237 821 447 844 662 150 7 2126.028 674 499 128 831 4 

400 4.499 751 495 177 639 267 396 02 14.998 829 822 856 441 699 608 131.237 821 447 844 662 151 2126.028 674 499 128 831 46 



TABLE II: The comparison of dipole sum rules, S^(n) and S^(n), for the H(ls) ground state. The 



exact expressions of sum rule are also presented in the fourth-column with 71 = \/ k 2 — a 2 Z 2 12]. 
The ratio ASi (n) / Sf xact (n) = [<S'f xact (n) — Sf (n)]/5f xact (n). The non-relativistic values are in the 
column 5^ R (n). The value of 5^(0) is not stable and gets smaller as the B-spline basis dimension 
is increased. The notation a[b] means a x 10 6 . 



Sum rule 



S+(n) 



5±(n) 



Sf*(n) 



Sf xact (n) 



ASi(n) 
5° xact (n) 



Si (-3) 
Si (-2) 

Si(-l) 
Si(0) 

%(1) 

■Si (2) 
Si (3) 



10.749 260 777 454 106 9 
4.499 751 495 886 496 76 

1.999 911 249 278 034 15 
0.999 955 631 350 807 45 



10.749 260 777 454 125 8 
4.499 751 495 177 639 27 

1.999 937 873 065 244 31 
l[-29] 



0.666 563 210 276 996 94 3.755 773 008 441 865 7y[4] 
1.298 802 722 313 -1.410 595 609 170 78 [9] 

5.298 0179 899 7[13] 



10.75 
4.50 

2.0 

1.0 

2 

3 
1 



(71 +1) (2 7 i + 1) 

3Z 2 




8 
3ofi 



2 /a 2 

4 

3a 4 \ 
2(7i 2 -l)(7i-2) 



71 + 



71(271-1) 



+ 3 



7[-19] 
-l[-29] 
2[-18] 

7[-16] 
-2[-12] 



influenced by the size of the confinement radius. The energy was converged to 25 significant 
digits for a basis with N = 400. 

Table [I] shows the convergence of the static multipole polarizabilities, af, for the H(ls) 
state as the dimension of the B-spline basis was increased from N = 100 to N = 400. The 
radius of confining cavity is R = 400 a.u. The static dipole polarizability af, is computed to 
a precision of 22 digits. The higher-order polarizabilities af, af, and a J have not achieved 
the same degree of precision, but are still computed to a precision of 21, 20, and 20 effective 
figures respectively. The present af = 4.499 751 495 177 639 267 396 02 a.u. is 4 x 10~ u 



a.u. larger than the result 4.499 751 495 142 92 a.u. of Goldman [13J. This is due to the 



different fine structure constant used. When the fine-structure constant a, is set to the value 
used by Goldman, namely 1/a = 137.035 999 074, the B-spline polarizability changed to 
af = 4.499 751 495 142 916 a.u. This is in perfect agreement with that of Goldman. All 
hydrogen atom sum- rules reported from now on use the N = 400, R = 400 a.u. B-spline 
basis. 

Exact expressions exist for a number of the dipole sum rules given by Eq. (|22|) . For 
example, the expressions for the exact non-relativistic electric-dipole sum rules Si(n) have 



been derived for n = —5, —4, • • • , 2 [30|,|33j. The non-relativistic dipole sum- rule diverges for 



n > 3. Expressions for some dipole sum rules for the Dirac hydrogen atom have been given 



by Drake and Goldman [12j. The Dirac equation sum rules were derived by using closure 
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TABLE III: The comparison of the H(ls) static multipole polarizabilities and sum rules with 



and without the negative energy states. Values for S, (0) are not numerically stable and tend to 
decrease as the basis is enlarged. The notation a[b] means a x 10 6 . 



Sum rule 






Non-relativistic 


S 2 (-3) 
<S 2 (-2) 

s a (-i) 
s 2 (o) 

Sa(l) 
5 2 (2) 
5 2 (3) 


26.747 582 450 922 508 1 
14.998 829 827 109 609 A 
8.999 384 961 848 033 62 
5.999 605 961 023 935 20 
4.799 574 122 244 155 50 
5.598 084 401 298 530 18 
20.043 653 259 626 


26.747 582 450 922 621 3 

i A nno onn ono occ a a i *7 

14.998 829 822 856 441 7 
8.999 544 703 293 683 54 
-l[-28] 

2.253 393 804 665 723 66 [5] 
-8.463 168 026 239 51 [9] 
3.178 616 122 279 928[14] 


26.750 

15.0 

9.0 

6.0 

4.80 

5.60 

20.80 


S"3 (-3) 
S3 (-2) 
Ss(-l) 
Ss(0) 

s 3 (i) 

S3 (2) 
S3 (3) 


204.041 400 069 326 002 
131.237 821 495 692 427 
89.992 366 277 948 754 1 
67.494 445 945 723 638 3 
57.851 751 717 231 297 2 
61.704 978 988 083 113 1 
100.225 824 655 056 308 


204.041 400 069 327 276 
131.237 821 447 844 662 
89.994 163 347 335 786 6 
6[-24] 

2.535 018 531 367 228 21 [6] 
-9.520 821 387 020 48[10] 
3.575 836 136 347 471 [15] 


204.06250 
131.250 
90.0 
67.50 

4jp = 57.8571428 
4|2 = 61.7142857 
^ = 100.2857142 


S 4 (-3) 
S 4 (-2) 

St(-l) 
S 4 (0) 

s 4 (i) 

S 4 (2) 
S 4 (3) 


3043.342 638 220 471 07 
2126.028 675 392 279 56 
1574.846 608 527 950 22 
1259.880 083 503 994 17 
1119.885 837 666 203 30 
1159.848 826 903 744 95 
1530.311 794 804 461 


3043.342 638 220 494 85 
2126.028 674 499 128 83 
1574.880 153 465 178 52 
9[-19] 

4.731 967 094 641 248 23[7] 
-1.777 189 196 560 64[12] 
6.674 763 648 090 144[16] 


3043.687 50 

2126.25 

1575.0 

1260.0 

1120.0 

1160.0 

= 1530.66666666 



to sum over the complete set of positive and negative energy states and the expressions are 
given in Table ILT1 The Dirac equation sum-rule for 51(3) is convergent due to cancellations 
between the terms with positive and negative energies. 

Table HT1 compares the dipole sum rules of the H(ls) with and without the contributions 
of the states in the negative energy sea. All the digits listed are converged with respect to 
further enlargement of the B-spline basis. The sum-rules, Sf~(0), S x (1) and S 1 (2), agree 
with the exact expressions to better than 15 digits. Agreement is not so good for Si (3) but 
in this case the sum is more sensitive to terms that occur at larger positive and negative 
energies. There was no evidence of convergence for Si (3) when the states of the negative 
energy sea were omitted from intermediate sum. This is consistent with the non-relativistic 
result of Lamm and Szabo (lij • 

The value S+(0), which omits the states from the Dirac sea, is close to the non-relativistic 
value of nuclear charge Z = 1. Upon making the substitution (p 2 ) = Z 2 in existing expres- 



sions 



we obtain the result 



5 +Levinger (0) = j _ + ^ ^ 

Evaluating this expression for Z = 1 gives, Si(0) = 0.999 955 6238, which is only 8 x 10 -9 
different from the B-spline evaluation. The degree of difference between Sf R (n) and (n) 
gets larger as n increases. The difference is 3 % for Si (2). 

The contribution that the negative-energy states make to the dipole sum rules depends 
on n. The negative-energy states of the Dirac sea contribute less than 2 x 10~ 5 to S^(— 1), 
2 x 10 -10 to S^(—2), and 2 x 10~ 15 to S^(— 3). This is not surprising. The negative energy 
states are located at energies of order —2c 2 . So the contributions of the negative energy 
states decrease as n in Eq. ( 122]) becomes increasingly negative. Conversely, the differences 
between the S x (n) and the S^(n) sum rules can be expected to increase as n increases. 
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TABLE IV: The static dipole polarizabilities for the ground-state of selected hydrogen-like ions. 
The present values are listed in the third- and fourth-columns for two sets of c = 1/a. All the 



tabulated digits of the present work are insensitive to further enlargement of the basis. The notation 
a [b] means a x 10 6 . 



z 


a+ 
"l 

r — 1 37 035 999 074 


at 

r — 1 37 035 999 074 








L37.035 989 5 




Goldman 


[13] 




1 


4.499 751 495 886 496 765 8 


4.499 751 495 177 639 267 4 


1 


499 


751 


495 


142 


915 


967 2 


4.499 751 495 142 92 


2 


0.281 187 875 627 153 384 5 


0.281 187 874 918 503 235 4 





281 


187 


874 


909 


822 


724 5 


0.281 


187 874 909 82 


5 


7.190 061 953 255 011 860[— 3J 


7.190 061 460 476 174 63[— 3] 


7 


190 


061 


244 


659 


087 


6 


-3] 


7.190 061 244 659 0[ 


-3] 


10 


4.475 171 382 242 160 041 [-4J 


iA7& 164 360 625 272 209[— 4] 


4 


475 


164 


357 


157 


090 


8 


-4] 


4.475 


164 


357 


157[- 


4] 


15 


8.778 661 031 860 895 31[— 5J 


8.778 591 625 838 392 08[— 5] 


8 


778 


591 


610 


447 


560 


3 


-5] 


8.778 


591 


610 


447[- 


5] 


20 


2.750 591 823 590 310 61[— 5J 


2.750 523 499 062 579 08[— 5] 


2 


750 


523 


490 


423 


618 


6 


-5] 


2.750 


523 


490 


424[- 


5] 


25 


1.112 456 189 324 034 04[— 5] 


1.112 389 181 457 920 41[— 5] 


1 


112 


389 


175 


944 


142 


1 


-5] 


1.112 


389 


175 


944[- 


5] 


30 


5.281 595 642 877 009 9[-6] 


5.280 940 730 404 758 7[-6] 


5 


280 


940 


692 


243 


592 


6 


-6] 


5.280 


940 


692 


243[- 


6] 


35 


2.798 031 223 308 353 9[-6] 


2.797 393 149 766 563 6[-6] 


2 


797 


393 


121 


842 


089 


7 


-6] 


2.797 


393 


121 


842[- 


6] 


40 


1.604 622 695 629 832 0[-6] 


1.604 002 839 548 263 7[-6] 


1 


604 


002 


818 


268 


128 


9 


-6] 


1.604 


002 


818 


268[- 


6] 


45 


9.767 839 136 814 269[-7] 


9.761 833 945 433 110[-7] 


9 


761 


833 


778 


188 


453 


1 


-7] 


9.761 


833 


778 


187[- 


7] 


50 


6.226 889 347 856 944[-7] 


6.221 086 480 106 640[-7] 


6 


221 


086 


345 


451 


685 


9 


-7] 


6.221 


086 


345 


451[- 


7] 


55 


4.116 918 654 470 464 [-7] 


4.111 325 157 474 914[-7] 


4 


111 


325 


046 


935 


820 


5 


-7] 


4.111 


325 


046 


936[- 


7] 


60 


2.802 469 149 798 750 [-7] 


2.797 090 474 417 353[-7] 


2 


797 


090 


382 


223 


3 13 


9 


-7] 


2.797 


090 


382 


224[- 


7] 


65 


1.953 091 120 155 380[-7] 


1.947 931 407 519 126[-7] 


1 


947 


931 


329 


604 


639 


5 


-7] 


1.947 


931 


329 


604[- 


7] 


70 


1.387 222 340 637 801 [-7] 


1.382 284 686 111 543[-7] 


1 


382 


284 


619 


529 


769 


8 


-7] 


1.382 


284 


619 


530[- 


7] 


75 


1.000 397 933 028 34[-7] 


9.956 846 315 732 27[-8] 


9 


956 


845 


741 


359 


051 


6 


-8] 


9.956 


845 


711 


359[- 


8] 


80 


7.301 102 574 925 93[-8] 


7.256 230 363 582 21[-8] 


7 


256 


229 


864 


059 


587 


8 


-8] 


7.256 


229 


864 


060[- 


8] 


85 


5.376 751 290 435 41 [-8] 


5.334 153 759 283 73[-8] 


5 


334 


153 


321 


793 


719 


5 


-8] 


5.334 


153 


321 


795[- 


8] 


90 


3.984 403 901 650 36[-8] 


3.944 093 881 570 48[-8] 


3 


944 


093 


496 


045 


404 


3 


-8] 


3.944 


093 


496 


045[- 


8] 


95 


2.962 871 397 452 loj-8] 


2.924 863 256 366 13[-8] 


2 


924 


862 


914 


773 


842 


2 


-8] 


2.924 


862 


914 


774[- 


8] 


100 


2.204 334 865 912 88[-8] 


2.168 647 587 493 68[-8] 


2 


168 


647 


283 


324 


507 


9 


-8] 


2.168 


647 


283 


325[- 


8] 



Table [TT1 shows that this indeed does happen. The difference between Sf(2) and Sf(2) is 
nine orders of magnitude. 

Table ILTT1 gives the sum rules for the higher-order multipoles for the hydrogen-atom ground 
state. The S^— 2), Ss(— 2), and S^— 2) are the multipole polarizabilities a2, 03, and 
respectively. The sum-rules, S^(n), omitting the states from the Dirac sea are within 0.1% 
of the non-relativistic values with the exception of S2 (3). This is also true for the sum-rules, 
Sf{n), with n < that also include the Dirac sea. 

The most striking results from Table [TTT1 are the Sf(0) sum-rules which do not exceed 
10~ 18 . Levinger et al 35| have pointed out that the Dirac Hamiltonian involves terms linear 
in the particle momentum p and that as a consequence the Bethe sum rule for exp(zq ■ r) 
should be identically zero. The expansion of exp(zq-r) implicitly involves dipole, quadrupole 
and octupole matrix elements. Therefore, it is expected that S^(0) = for all i. 

The contributions of the negative-energy Dirac Sea to the Sf{n) sum-rules are actu- 
ally greater than the contributions from the physical states for n > 1. They exceed the 
contribution from the physical states by amounts from 4 to 14 orders of magnitude. 



B. Polarizabilities for the hydrogen isoelectronic series 

Table IIVI presents the static dipole polarizabilities for a number for hydrogen-like ions 
in their ground-state. All the digits listed in this table are converged with respect to fur- 
ther enlargement of the B-spline basis. In order to facilitate comparison of the present 
polarizabilities with those of Goldman (l3j|, we repeated the calculations but used the the 
same speed of light c = 137.035 989 5 (in atomic units) as Goldman. The agreement 



9 



with the polarizabilities of Goldman could hardly have been better. At Z — 10 we got 
af = 4.475 164 357 157 0908 x 10~ 4 a.u., in agreement with all published digits of Goldman. 
The same is true for the polarizability at Z = 90, namely 3.944 093 496 045 4043 x 10~ 8 
a.u. This level of agreement was achieved for all values of Z from 1 to 100. The only 
disagreements amounted to ±1 in the last significant digit reported by Goldman (l3| . 



0.0 
x -0.2 

a 

5 -0.4 

I 

+1 — I 

3 -0.6 



20 40 60 80 100 
Nuclear charge Z 




FIG. 2: (color online). The impact of relativistic effects on the multipole polarizabilities for the 
hydrogen isolectronic series. The ratio (af — af R )/a^ R is plotted. 

The higher-order polarizabilities, af of the ground-states of some selected hydrogen-like 
ions are presented in Table |Vj All the reported digits are insensitive to further enlargement 
in the B-spline basis. Fig. [2] shows the influence of relativistic effect on multipole polariz- 
abilities. The relativistic effect becomes larger as the nuclear charge, Z, is increased. The 
relative size of the relativistic effect is smallest for the dipole polarizability and largest for 
a^. 




FIG. 3: (color online). Plot of 10 9 Aai/Z 2 as a function of nuclear charge, Z. 

The difference of the af and af polarizabilities from the Kaneko polarizabilities are 
illustrated in Fig. [3j We define Aaf = (af — af) with a similar relation used to define 
Aaf. Fig. [3] plots 10 9 Aai as a function of Z. These are seen to go to a constant value as 
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Z — > 0. From Eq. ( 120]) we deduce 



Aaf 



2Z A 



A 2 + |)(nZ)- • Ofn'Z'^ 



(24) 

§. Fig. m 



This expression can only go to a constant in the Z — > limit when A2 = 
demonstrates that a^, aj 1 ", and af are equal to order 0(a 2 Z 2 ). The different Z — > 
asymptotes for Aa^ and Aa| indicate that the 0(a 4 Z 4 ) terms are different for and af. 
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FIG. 4: (color online). Plot of 10 



15 'af 



)/Z 2 as a function of nuclear charge, Z. 



Expressions for af giving terms to 0(a 4 Z 4 ) J9j, 0(a 6 Z 6 ) [10| and to all orders 11] have 
been derived. The 0(a 4 Z 4 ) expression of Zon [9] is 



a-, 



Zon 



9 



2Z A 



, 28 2 31 + 27T 2 

1 ttZ H 

27 V 7 432 



« 4 Z 4 ) + 



(25) 



Figure H] plots 10 15 (a 1 z — af on )/Z 2 as a function of Z. The Z — )■ limit of this difference 
demonstrates that the present calculations are in agreement with the analytic expression to 
0(a 4 Z 4 ). This provides a mutual validation of the B-spline calculations and the analytic 
expressions. 

Figure [5] plots 10 9 Z 4 Aa 3 as a function of Z. These are seen to go to a constant value 
as Z — y 0. By an analysis similar to that performed for the dipole polarizability, one can 
deduce that af, af, and af are equal to order 0(a 2 Z 2 ). The different Z — > asymptotes 
for Aaf and Aaf indicate that 0{a 4 Z A ) terms are different for af and af. 

The analysis demonstrating that the differences between af and af only appear at 
0(a 4 Z 4 ) has already been reported It has previously been shown that these polar- 
izabilities are agreement with the Kaneko 0(a 2 Z 2 ) It is also possible to plot 10 9 Z 6 Aa4 
as a function of Z giving plots similar to Figs. [3] and [5j This demonstrates that af and af 
agree with af at the 0(a 2 Z 2 ) level and the difference between af and af occurs at the 
0(a 4 Z 4 ) order. 



C. Sum-rules for the hydrogen isoelectronic series 



The non-relativistic TRK sum-rule, Si(0) gives a value of unity [30|, |33|, |37( for all 
hydrogen-like atoms and ions. However, ^(O) is exactly zero while the sum-rule, Sf(0) 



11 



TABLE V: Relativistic multipole polarizabilities (in a.u.) for the ground states of the hydrogen 
isoelectronic series. All the figures listed are accurate. The notation a[b] means a x 10 b . 



7 




a 3 




1 


14.998 829 822 856 441 699 


t O t £ ~\ o T C\ C\ 1 A A T C~\ A A /~* /~* C\ 

131.237 821 447 844 662 


2126.028 674 499 128 83 


2 


0.234 301 867 935 791 210 


0.512 505 037 523 770 47 


2.075 551 5460 612 051 9 


5 


9.581 285 372 324 0453 92[— 4J 


o oro O 1 ( \ poo TOT A1 z 1 oT A~\ 

3.352 210 608 787 016 2 [ — 4J 


2.171 618 426 945 541 1[— 4J 


10 


1.488 319 383 913 4ll 04[— 5J 


1 o oo o r* o ooo tot o o ,i T f~> I 

1.300 352 899 787 624[— 6J 


O -| O ,1 -|OT /"* ,1 r" TKO Ol Tl 

2.104 187 645 750 314[ — 7J 


15 


1.293 852 351 688 892 4[— 6J 


r* /~i -i a o tt /loo oot otT nl 

5.014 877 480 967 07|— 8J 


o /"i o 1 r*oo or*^ r*o-i r*T ol 

3.601 503 954 501 5[— 9J 


20 


o rtw 1 i o r" o o ata to o T Tl 

2.271 146 583 050 793[— 7J 


A ooo o A o oto o y» o oT ol 

4.938 640 072 269 2[— 9J 


-\ OOI 0<^*0 A A O 01TT iol 

1.991 062 443 017| — 10J 


25 


r* o T o^r' r* o r* tot o o T nl 

5.847 845 585 737 33[— 8J 


o -i -i o o r* o ~i / * iT a oooT i aI 

8. HO 859 162 392[-lOJ 


O O O T O TO TT"I OOT -1 -1 1 

2.087 370 771 99[— llj 


30 


■t o -i r" f 1 k to -i o o r rnf ol 

1.915 515 761 865 58[— 8J 


-i oot oo/ 1 /too /irnf i aI 

1.837 296 630 650[— 10J 


O OTO ~1 OO P* O 1 tT -i ol 

3.273 123 521 7[— 12J 


35 


TOOT /(TO O ^ r r n A i [ ol 

7.397 473 245 589 1|— 9J 


r - io/ 1 oto oto r* r\\ i 1 1 

5.186 973 978 69[— llj 


r> Tr*o lor - r/'Af 10I 

6.763 105 560[— 13J 


i o 

40 


O C\ 't O ki*^/- OTO ri/i A O T O 1 

3.218 326 876 369 0[— 9j 


1 T1 T /^Tl -t -t f* TO T 111 

1.717 671 116 72|— llj 


1 TOT O /"> T OOTr 1 ol 

1.707 067 337[— 13J 


45 


1 roi " ( ' \ v r\c\ 010 T T nl 

1.531 561 509 916 7[—9\ 


f* /ii r oo A o A o i r i ol 

6.415 324 043 1 [ — 12J 


pf 1 1 ooo oof "1 /ll 

5.011 809 33[ — 14J 


r o 

50 


7.812 859 401 235[-10J 


O f O O i^OO f T 1 oT "1 Ol 

2.630 602 571 9[— 12J 


-| f A O O 1 OTP -f Al 

1.654 931 37[— 14J 


55 


4.210 472 655 409[— 10J 


1.161 555 467 5[ — 12J 


5.999 556 2[— 15J 


60 


O O T -| 1 J »7 Afl A J i( T "I ol 

2.371 147 053 044[— 10J 


r" ii /I fi r'TO oooT i ol 

5.443 579 080|— 13J 


o o i r ooo oT 1 P*l 

2.345 208 2 [ — 15J 


65 


i ooo /"'it rrr <i of i ol 

1.383 617 655 412[— 10J 


O/^TT /IP'T /lOOT 1 ol 

2.677 457 400[— 13J 


o T/io Aor f ~\ r*~\ 

9.748 095[— 16J 


70 


8.309 087 512 23[-ll] 


1.369 821 733[-13] 


4.261 037[-16] 


75 


5.106 469 950 92[-ll] 


7.235 969 19[-14] 


1.940 914[-16] 


80 


3.196 013 748 39[-ll] 


3.921 694 89[-14] 


9.141 67[-17] 


85 


2.028 253 121 49[-ll] 


2.168 463 36[-14] 


4.421 83[-17] 


90 


1.299 794 490 85[-ll] 


1.216 900 77[-14] 


2.182 71 [-17] 


95 


8.376 878 675 0[-12] 


6.895 117 0[-15] 


1.092 81[-17] 


100 


5.405 559 183 5[-12] 


3.923 335 2[-15] 


5.514 2[-18] 



is almost equal to 1. The appropriate method to choose for the evaluation of the TRK 



sum rule has generated considerable discussion |l2|, |35|, |36|, |38f|43| . Table [VI] compares the 
present B-spline values of (0) and compares them against the earlier calculation of Drake 
and Goldman [l2j]. Keeping in mind the limited precision of the earlier calculation, the 
agreement with the Drake and Goldman calculation is perfect. 
Figure E] shows 

ASi(0) = S+(0) - S+ Levinger (0), (26) 

plotted as a function of Z. It is noticed that ASi(0)/Z 4 goes to a constant as Z — > 0. 
This demonstrates that the present 5*^(0) is in agreement with the expression of Levinger to 
order 0(a 2 Z 2 ). It also demonstrates that the next term in the expression for (0) occurs 
at the 0(a 4 Z 4 ) level. The near linear behaviour of AS'i(0)/Z 4 at small Z indicates that the 
expansion for (0) contains a term of 0(a 5 Z 5 ^ 



While the B-spline calculations of 5* + (0) are compatible with 0(a 2 Z 2 ) expressions |35l 136 
they cannot be reconciled with the 0(a 4: Z 4: ) expression of Cohen. A simple analysis near 
Z = suggests that 

S+(0) w 1 - + 2.71 a 4 Z 4 - 6a 5 Z 5 + 0(a G Z 6 ) . (27) 
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TABLE VI: Comparison of the S± (0) sum rules. All digits 


are stable with respect to further 


enlargement of the B-spline 


basis. 




Z 


Present 


Ref. [12] 


1 


0.999 955 631 350 807 


0.999 9556 


2 


0.999 822 612 102 297 




5 


0.998 894 823 187 627 




10 


0.995 622 481 263 678 


0.995 62 


15 


0.990 287 581 618 103 




20 


0.983 023 671 163 131 


0.9830 


25 


0.973 973 703 862 452 




30 


0.963 278 628 607 378 


0.9633 


35 


0.951 070 787 251 835 




40 


0.937 470 188 595 043 


0.9375 


45 


0.922 582 481 520 977 




50 


0.906 497 887 620 449 


0.9065 


60 


0.871 018 387 592 671 


0.8710 


70 


0.831 424 017 561 149 


0.8314 


80 


0.787 815 483 542 815 


0.7878 


90 


0.739 933 345 752 064 


0.7399 


100 


0.686 987 401 548 771 


0.69 



It has not been possible to reconcile the coefficient of 2.71 with Eq. (8) of Cohen |36j. But 
it is unclear how to interpret (p 4 ) of Eq. (8) in j36|. The plot of 5*^(0) depicted in Fi g. 3 
of 36[ is certainly compatible with the present B-spline calculation. However, Fig. 3 of 36l 
plots the 0((p 2 )) approximation to 5^(0) and this is certainly not equal to 1 — 5a 6 z . 
Figure [7] shows the difference 

ASi(n) = St(n) - S+{n) , (28) 
plotted against Z for n = —1 and n = —3. It is noticed that limz_> ASi(— 1) and ASi(— 3) 
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FIG. 6: (color online). Plots of 10 9 AS\(0) / 'Z^ versus nuclear charge, Z. 
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FIG. 7: (color online). Plots of AS\{— 1) and A5i(— 3) versus nuclear charge, Z. 

both go to a constant as Z — > 0. Figure |3] established that ASi(— 2) also has the same 
Z — > limiting behaviour. Writing either of the S\(n) in the form 

S 1 (n) = Sf R (n) (1 + c 2 a 2 Z 2 + c 4 a 4 Z 4 + . . .) , (29) 

allows one to deduce that the c 2 coefficients are different for Sf~(— 1) and Si(— 1) since 
5*^ R (— 1) = 2/Z 2 . However, one deduces that the c 2 and C4 coefficients are actually the same 
for Sf~(-3) and S?(-3) since S? R (-3) = 43/(4Z 6 ). 

D. Analytic expressions for the multipole polarizabilities of hydrogen-like ions 

Analytic expressions were derived for af by performing a least squares fit of the polariz- 
abilities to Eq. fl20|) . The polarizabilities were divided by the non-relativistic values prior to 
the fit. The value of A2 was fixed at the values of Kaneko for £ = 1,2 and 3. The A2 value 
for i = 4 was determined by evaluating Eq. (36) of Ref. p}. The value of A 4 for af was 
set to the value from Zon [9j . Table IVHI lists the numerical values of X 2 i coming from the 
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FIG. 



function of nuclear charge, Z. 



(color online). Relative difference between a« and the fits to this using Eq. (|20p as a 



TABLE VII: The Eq. (|20p fits to the multipole polarizabilities of hydrogen isolectronic series ground 
states. 



Term 



a 



± 



[13] 



n i 



a 



a.. 



a NR 


9 


9 


15 


525 


8505 




2Z 4 


2Z 4 




4Z 8 


4Z 10 


A 2 


28 


28 


879 


5123 


33251 


~27 


~27 


~600 


2940 


17010 


A 4 


0.117 451 87(1) 


0.117 451 870 668 402 


0.502 471 315 


0.854 144 263 


1.177 235 432 


A 6 


0.007 482(1) 


0.007 692 784 


-0.014 151 521 


-0.102 874 518 


-0.228 232 960 


A 8 


0.0010(1) 


-0.003 271 333 


0.002 052 103 


0.001 434 636 


0.011 629 938 


Aio 




0.006 117 861 


-0.000 261 805 


0.001 019 239 


-0.000 189 306 


Al2 




-0.013 528 604 









a 



a i 



A 2 

A 4 
A 6 
A 8 
Aio 



_28 
~27 
0.171 953 291 
-0.069 671 936 
0.075 248 612 
-0.051 443 668 



879 

~600 

0.601 241 304 
-0.164 544 588 

0.118 942 451 
-0.057 448 561 



5123 

~2940 
0.981 404 521 
-0.321 798 414 
0.171 613 014 
0.073 735 494 



33251 
~ 17010 
1.323 923 421 
-0.506 146 068 
0.236 007 349 
-0.093 122 365 



fit. These coefficients give a more precise representation of the exact dipole polarizabilities 
than two previous representations 12, 1^. The expressions for the quadrupole and octupole 
polarizabilities are novel. 

The quality of the fit to the B-spline af can be seen from Fig. El The quality of the fits 
are of very high accuracy at the smaller values of Z. This occurs since the leading X 2 a 2 Z 2 
term uses the exact value of A 2 . The quality of the fit is degraded at larger values of Z. 
However, the maximum relative error in the analytic expressions only exceeds one part per 
million for values of Z close to 100. 

Equation ( 120]) was also used to create an analytic expression for a[ . In this case, the \^% 
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TABLE VIII: The second-order dispersion coefficients (in a.u.) for the H(ls)-H(ls) and H(ls)- 
He + (ls) systems. Results are given for the sum rules evaluated with and without the states of the 
Dirac sea. All tabulated digits are accurate. The notation a [b] means a x 10 b . 









Non-relativistic 






H(ls)-H(ls) 




c 6 


6.498 392 250 007 09 


6.498 392 245 754 06 


6.499 026 705 405 84 


C 8 


1.243 840 307 694 35[2] 


1.243 840 306 577 93 [2] 


1.243 990 835 836 22[2] 


Cio 


3.285 370 791 861 60[3] 


3.285 370 788 289 10[3] 


3.285 828 414 967 42[3] 






H(ls)-He+(ls) 




c 6 


0.657 548 755 759 311 


0.657 548 758 416 787 


0.657 716 656 238 770 


C 8 


8.335 406 342 724 22 


8.335 406 384 081 09 


8.337 819 589 166 31 


Cio 


1.588 773 716 759 79 [2] 


1.588 773 725 479 53 [2] 


1.589 267 575 526 71 [2] 



parameters with i > 1 were treated as fitting parameters. The results of the fit are tabulated 
in Table ECD 



IV. DISPERSION COEFFICIENTS 

The long-range dispersion interaction between two spherically symmetric atoms can be 
written 

oo „ 

W*)~-£^' (30) 

n=3 

The dispersion coefficients, C*2 n can be evaluated using oscillator strength sum rules. The 
explicit expression is 

r ( 2 -- 2 ) ! 



=1 L=l 



4(2Z i )!(2/ i ) 



x 



Ati) Ah) 

T J Y Jn -" -. (31) 



where £i + £j + 1 = n and £A,gi is the excitation energy from state g to state i for atom A. 
The sum implicitly includes the continuum, and j^Ag\ ^ s ^ ne oscillator strength of multipole 
£i connecting the state g to the excited state % for atom A. Considerations of molecular sym- 
metry do not have a direct effect on Eq. (13"T|) when both atoms are in spherically symmetric 
states. 

It is surprising that there has not yet been any calculation of the hydrogen dimer disper- 
sion coefficients based on oscillator strengths from the Dirac equation. This is rectified in 
Table IVLTfl where the C$, Cg and Cio coefficients are given for two hydrogen atoms in their 
ground states. Table IVIIII also gives the dispersion coefficients between a hydrogen atom 
and a He + ion. 

The use of the Dirac equation leads to the H-H C 6 being reduced by 0.00063 a.u. or 
0.0098%. The relative difference is about twice as large as the difference between the rela- 
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tivistic and non-relativistic polarizabilities. The reduction in the size of C% is larger for the 
H-He + system, being about 0.026%. 



V. CONCLUSIONS 

A computational investigation based on B-spline methods has been used to investigate 
the polarizabilities and related sum rules of the hydrogen isoelectronic series. Dipole polar- 
izabilities have been computed to a higher precision than any previous calculations. One 
distinction with previous calculations is that results were also reported for calculations where 
the negative-energy Dirac sea is excluded from the intermediate sum. The agreement with 
previously derived analytic expressions 0, 0, 10] for the dipole polarizability could not be 



better. High precision calculations of the multipole polarizabilities for i = 2,3,4 are also 
given. The present results provided a computational validation of the earlier works of Kaneko 
[7J and Zon [9f. The af polarizabilities are in agreement with the Kaneko expressions at 
the 0(a 2 Z 2 ) level. The af polarizability is also in agreement with the expressions of Zon 
]9]. which includes terms at the 0(a A Z A ) level. 

Precise values for other oscillator strength sum rules have also been computed. The 
sum-rule 5^(0) = Ylifgf = provides a valuable check of the numerical reliability of the 
calculations. The sum- rule, <SV~(0) has been shown to be compatible with the 0(a 2 Z 2 ) 
expression of Levinger et al [35]. It is also compatible with earlier numerical calculations 



12, 36 



One aspect of the present work that represents a departure from earlier work has been the 
treatment of the states of the negative energy sea. Existing practice is that calculations of 
polarizabilities include the states of the negative energy sea, while calculations of the Bethe 
sum rule tend to omit these state. The philosophy of the present work has simply been to 
do two calculations for most properties, those that include the states of the Dirac Sea and 
those that omit them. 

Analytic expressions for af and ol\ to relative precisions not exceeding 10 -6 have been 
obtained by fitting an (aZ) n expansion to the computed polarizabilities. The Cq, C$ and C\q 
dispersion coefficients for the long-range H-H and H-He + interactions were also computed. 
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